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On  the  Lagrangian  Description  of  Vorticity 

J.  Casey  &  P.  M.  Naghdi 

Communicated  by  D.  Owen 


Abstract 

This  paper  exploits  a  Lagrangian  form  of  the  vorticity  field  that  is  based  on 
a  formula  of  Beltrimi,  the  merits  of  which  do  not  seem  to  have  been  fully  appre¬ 
ciated.  A  number  if  new  results  are  readily  derived  by  using  the  Lagrangian  de¬ 
scription,  and  the  classical  vorticity  theorems  are  also  included.  The  results,  being 
purely  kinematical,  apply  to  all  deformable  media  (including  viscous  fluids)  and 
may  also  be  of  value  in  computations. 

1.  Iatrodactfoa 

Although  it  has  been  rarely  discussed,  the  Lagrangian  description  of  vorticity 
does  offer  some  advantages  over  the  conventional  (Eulerian)  description.  Thus,  for 
example.  Truesdell1  has  shown  'hat  the  division  of  vorticity  transport  into  con¬ 
vective  and  diffusive  mechanisms  becomes  especially  transparent  when  a  Lagran¬ 
gian  viewpoint  is  adopted.  Furthermore,  as  will  become  apparent  below,  the  treat¬ 
ment  of  conditions  for  the  materiality  of  vortex-lines  in  a  deforming  continuum 
is  particularly  straightforward  when  phrased  in  Lagrangian  terms.  Because  of 
these  and  other  reasons  of  a  more  general  nature,  it  seems  to  us  that  a  Lagrangian 
description  of  vorticity  should  be  more  fully  investigated  than  it  has  been. 

Since  the  history  of  the  subject  is  quite  tangled,  it  is  desirable  to  provide  some 
background  information.  In  the  context  of  classical  hydrodynamics,  Cauchy 
showed  that  under  certain  circumstances  the  current  vorticity  vector  at  at  a  fluid 
particle  could  be  explicitly  written  in  terms  of  its  initial  value1  too0.  One  consequence 

1  See  sect.  84  of  Truesdell’s  monograph  {1]  and  also  an  earlier  paper  by  him  [21_ 

1  See  our  eqn.  (3.25),  in  which  F  signifies  the  deformation  gradient  tensor  and  J 
is  its  determinant.  Cauchy’s  formula,  presented  in  1813,  appears  as  eqns.  (16)  on  p.  40 
of  [3],  As  Truisdsll  has  pointed  out  in  his  discussion  [1,  p.  173ff.],  the  significance 
of  Cauchy's  formula  was  not  appreciated  in  the  decades  following  its  publication.  A 
notable  exception  is  Stokes,  who  became  increasingly  aware  of  its  importance  (see  [4, 
p.  108]  and  [5]).  A  discussion  of  Cauchy’s  formula  is  given  in  Lams  (6,  p.  205]  and  it  is 
treated  also  in  Serrin  (7  p.  1421.  Truesdell  &  Tourin  [8,  p.  421],  and  Batchelor 
(9.  p.  2761. 
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of  Cauchy’s  formula  is  that  vortex-lines  in  a  circulation-preserving  flow  are  carried 
as  material  lines.  In  contrast,  in  a  more  general  motion,  the  vortex-lines  no  longer 
move  as  material  lines,  and  Cauchy's  formula  no  longer  holds.  Nevertheless, 
it  is  still  possible  to  have  an  expression  of  the  same  form  as  Cauchy’s  in  order  to 
associate  with  the  current  vorticity  field  a  time-dependent  field  a* 3  defined  on  the 
reference  configuration  of  the  continuum.  To  our  knowledge,  Beltrami  [10] 
was  the  first  to  give  such  a  formula,  and  it  appears  to  have  gone  unnoticed  ever 
since.3 5 * *  Beltrami's  development  is  rather  convoluted,  but  he  definitely  has  a  La- 
grangian  vorticity  vectoi  ana  ne  also  deduces  a  few  of  its  properties.  In  the  La- 
grangian  representation  of  vorticity  which  Truesdell  presents  in  [1],  Beltrami's 
vorticity  vector  is  not  explicitly  introduced,  but  a  large  number  of  results  are 
derived  in  Lagrangian  form.*  In  the  present  paper,  we  employ  Beltrami’s  formula 
as  the  cornerstone  for  the  construction  of  the  vorticity  field  in  Lagrangian  form. 
The  discussion  is  purely  kinematical  and  >  herefore  applies  to  all  deformable  media 
including,  for  example,  viscoelastic  and  elastic-plastic  materials,  as  well  as  viscous 
fluids. 

The  contents  of  the  paper  are  as  follows.  In  Section  2,  we  briefly  review  the 
Eulerian  description  of  the  vorticity  field.  It  should  be  noted  that  in  the  Eulerian 
description  the  curl  of  the  acceleration  is  a  relatively  complicated  function  of 
the  vorticity  in  that  it  also  involves  the  rate  of  deformation  tensor  (see  (2.12) 
and  (2.9)i  :).  The  Lagrangian  representation,  which  is  introduced  in  Section  3, 
furnishes  a  relationship  between  acceleration  and  the  vorticity  that  could  not  be 
simpler:  the  curl  of  the  acceleration  is  the  material  derivative  of  the  (time- 
dependent)  Beltrami  vorticity  vector  (see  (3. 19)t  and  (3.17,)).  Many  theoretical 
results  follow  directly  from  this,  and  it  may  also  be  of  computational  value.  In 
Section  4.  Eulerian  and  Lagrangian  representations  for  circulation  are  discussed, 
and  the  time-evolution  of  vortex-tubes  is  analyzed.  In  Section  5,  we  present  ne¬ 
cessary  and  sufficient  conditions  for  vortex-lines  to  be  transported  as  material 
lines.  Also  included  in  Sections  4  and  5  are  the  connections  between  the  present 
results  and  classical  ones,  especially  the  vorticity  theorems  of  Helmholtz  [11] 
and  Kelvin  [12].s 


2.  Preliminaries.  Eakrian  description 

Consider  a  deformable  continuum  9  moving  in  three-dimensional  space  and 
let  If  be  a  typical  particle  of  3t.  Also,  let  X  and  x,  respectively,  denote  the  position 
occupied  by  X  in  a  fixed  reference  configuration  *°  of  &  and  in  the  present 
configuration  »  at  time  t.  The  motion  of  &  is  described  by  the  mapping  x  =  %(X,  t). 
The  particle  velocity  v,  deformation  gradient  F,  determinant  of  F,  and  right 


3  The  formula  is  represented  by  eqns.  (13)  of  (10}.  These  correspond  to  eqn.  (3.2), 
below,  when  written  in  its  equivalent  form  /<*  »  Fat9. 

*  Although  Truesdell  cites  [10]  in  sect.  84  of  [1],  it  seems  that  he  overlooked  the 
value  of  Beltrami’s  eqns.  (13). 

5  Various  statements  and  proofs  of  these  theorems  can  be  found  in  [1,  6,  7,  8,  9], 

and  the  reader  may  also  wish  to  consult  (13,  Sect.  11.13]  and  [14,  Sect.  11].  A  discussion 

of  the  continuing  usefulness  of  the  classical  results  can  be  found  in  [15]. 
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Cauchy-Green  measure  of  deformation  are 

/  =  det  F  >  0,  C  =  FtF, 


d%  r.  d X 

v  =  x  —  —  F  =  — 

dt  ’  3JT’ 


(2.1) 


while  the  velocity  gradient,  rate  of  deformation  tensor,  and  vorticity  tensor  are) 
given  by 


cv 

£=-,  D  —  Liym  =  $  (L  +  L1),  W  =  L,kw  =  i(L  —  L7).  (2.2) 

We  recall  that 

F  =  LF,  ~F~l  =  -F~lL,  J  =  /trD  =  /divv.  (2.3) 

The  vorticity  vector  at  and  tensor  W  satisfy  the  relations 

at  =  —e[W]  =  curl  v,  W  =  —  \  set,  2Wc  =  a>  x  c,  (2.4a) 

where  e  is  the  permutation  tensor  and  c  is  an  arbitrary  vector.  Relative  to  a  fixed 
orthonormal  basis  {«,;  i  =  1,  2,  3},  (2.4 a), i2,j  have  the  representations 


a)  i - e,jk  WJk  —  £:jkvkj,  Wtj - isiJka)k, 


(2.4  b) 


w  here  the  usual  convention  of  summation  over  repeated  indices  has  been  employed, 

c() 

and  where  ( )j  =  — — .  We  observe  that 


dxj 

Wat  —  0,  Lai  =  Dai  —  LTai,  divco  =  0. 


(2.5) 


A  motion  is  irrotational  at  time  t  if  and  only  if  at,  or  equivalently  W,  vanishes  at 
time  t  for  each  particle  X  of 

The  acceleration  and  its  spatial  gradient  are 


and  it  is  clear  that 

L  =  A  —  L2  =  A  -  D2  —  W2  -  DW  —  WD, 

D  =  —  D2  —  W2,  W  —  A,kw  —  DW  —  WD, 


(2.6) 


(2.7) 


where  the  symmetr  c  and  skew  parts  of  A  are  defined  by  formuae  of  the  type  (2.2)3 

A 

Recalling  the  formula*  for  the  convected  rate  W  of  W  and  noting  that  WTW  = 
—  W2,  we  have 


W  —  W  A-  LTW  +  WL  =  A,kw. 


(2.8) 


‘  See  Truisdell  Sl  Noll  [16,  p.  97].  Note  that  a  tensor  does  not  have  to  be  objective 
for  a  formula  of  the  type  (2.8)  to  apply. 
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Just  as  with  the  skew  tensor  W,  we  may  uniquely  associate  a  vector  %  with  the  skew 
part  of  the  acceleration  gradient.  Thus, 

*  =  -e[A,kw]  =  curia,  A,kw  =  —\e\.  (2.9) 

Expressing  v  as  a  function  of  the  spatial  variable  s  and  r,  one  readily  obtains 

tv 

a  —  —  -  rox  v  -t-  i  grad  (v  •  v),  (2.10) 

and  hence 

cut  tut 

'  =  — — f-  curl  (utxv)  ~  — — 1-  div  (ut  ®  v  —  v  to),  (2. II) 
Ct  of 

where  ®  denotes  the  tensor  product  of  vectors7  and  div  (cu  g  v)  =  et. 

Consequently,  in  view  of  (2.5)2,4. 

*  =  ti*  ■+■  to  div  v  —  Deo.  (2.12) 

If,  for  an  arbitrary  vector  u,  one  defines  the  Tnesdell  rate 8 9  of  u  by 

u  =  ti  +  {(tr  D) /- L}u,  (2.13) 

it  then  becomes  dear  that 

*  —  at.  (2. 14) 


3.  Ligrufiia  description 

Recall  that  the  adjugate  F *  of  F  satisfies  the  relations 

F*  =  /(F~ ‘)r,  (Fa)  x  (Fft)  =  F*(a  x  6),  det  F*  =  /J  >  0.  (3.1) 

for  any  vectors  a,  b.  We  define  time-dependent  vectors  to3  and  \°  by* 

to3  =  JF-'ut  =  jt-‘  (to),  (3.2) 

=  *~l  {*}•  (3.3) 

As  will  quickly  become  evident  from  its  properties,  the  vector  to3  provides  a 
Lagrangian  description  of  vorticity.  In  component  form,  the  relationships  between 

7  Recall  that  (a  ®  b)  c  —  ab  ■  e  for  any  vectors  a,  b,  c.  t 

*  This  terminology  is  introduced  to  reflect  certain  properties  which  the  rate  u  has 
in  common  with  the  Truesdell  rate  of  an  objective  tensor  (see  Casey  Sl  Nachdi  [17, 
p.  357]). 

9  The  operation  .f'  is  known  as  a  Piola  transformation.  Operations  of  this  type  com¬ 
monly  appear  in  transformation  formulae  connecting  Lagrangian  and  Eulerian  descrip¬ 
tions  of  fields  (e.g.,  the  stress  field).  The  vorticity  vector  «*°  (actually,  was  intro¬ 
duced  by  Beltrami  [10],  who  used  a  procedure  somewhat  different  from  that  given 
here. 
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co3  and  to  are 


ojAa  =  JXA<i<jiit  a)'  =  —xIA0)A°, 


(3.4) 


where  coAa  are  the  components  of  tu°  relative  to  a  fixed  orthonormal  basis  EA 
(which  may  or  may  not  coincide  with  et).  Similar  equations  hold  between  the 
components  <x,  and  xA°. 

Observing  the  identity, 


=  0,  or  Div  F*  =  0, 

(3.5) 

we  deduce  from  (3.4),  that 

Div  o»°  =  a)ArA°  =  JXAto)l'A  =  J  div  co. 

(3.6) 

Hence,  in  view  of  (2.5)*, 

Div  co°  =  0. 

(3-7) 

Letting 

W3  =  -*MO°, 

(3.8) 

we  note  that  formulae  of  the  type  (2.4a),  *  and  (2.5),  also  hold  for  fF°andto°. 
For  any  vector  c,  we  have 


to  x  c  =  (y  Fco° )  x  F(F"‘  c)  =  y  F*(w3  x  (F-1  c)) ,  (3.9) 

where  use  has  been  made  of  (3.2),  and  (3.1)2.  Then,  applying  (2.4 a)4  and  its  ana¬ 
logue  for  the  Lagrangian  quantities,  and  invoking  (3.1),,  we  find  that 

W°  =  FTWF.  (3.10) 


We  also  observe  that  irrotational  motions  are  characterized  by  the  vanishing  of 
to 3  (or  equivalently,  W°). 

In  anticipation  of  later  results,  we  transform  the  velocity  and  acceleration 
fields  by  means  of  the  formulae10 


and  also  let 


=  Ftv, 

a°  =  FTa, 

(3.U) 

8v° 

da° 

(3.12) 

~Tx' 

A°  _  - 

A  8X ' 

It  can  then  be  easily  verified  that 

L°  =  FTLF  +•  xtyAgotEA  ®  Eg, 

A0  =  FTAF  -J-  xl<AgfllEA  ®  Eg, 

i>A°  =  +  aA°. 

LAb°  —  (^°)j  —  puiei  +  +  AAg  . 


10  It  is  worth  pointing  out  that  the  transformation  (3. 11),  could  have  been  introduced 
prior  to  (3.2),.  The  connection  between  these  two  fields  will  appear  in  (3.15)*, 
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Breaking  L°  and  A°  in  (3.13),  2  into  their  symmetric  and  skew-symmetric  parts, 
we  obtain 


Lsym°  =  FrDF  J-  xlABv,EA  ®  £s,  L,kw°  =  FTWF  =  r, 

■4,rm 1  —  FTAiymF  4-  x,ABalEA  ®  Eg,  A,iw°  =  FrA,kmF , 


(3.14) 


where  use  has  een  made  of  (2.2)j  4  and  equations  of  the  type  (2.2)j,j  and  of  (3. 10). 

It  follows  from  (3.11),,  (2.2),i2i4,  (3.14)3,  and  a  formula  of  the  type  (2.4b), 
that 


Curl  v°  —  £ukvk ,/Ei  —  S/jgXf  gXj  jW^Ef 

—  B.jgWg/E, 


(3.15) 


Hence,  in  view  of  (2.4a)2  and  (3.2)  2, 

Curl  v°  =  .«*-»  (curl  v}. 

Likewise, 

Curl  o°  =  s°  =  {curl  a}. 

Next,  we  deduce  from  (3.13)3  that 

eUk(.vK°).J  —  etJKvlJUVl  +  elJKvU°t,K  +  EfJKaK.J° 
—  s!JKaK./ » 


(3.16) 


(3.17) 


(3.18  a) 


Curl  v°  =  Curl  a°. 

It  follows  from  (3.15)4,  (3.18b),  (3.17),,  (3.3),  and  (2.14)  that 

=  *°  =  *-l{«}.  =  Aikw°, 


(3.18b) 


(3.19) 


where  a  formula  of  the  type  (2.9)3  has  been  used  in  the  last  step. 

Suppose  that  at  some  time  t  =  f0,  the  body  occupies  its  reference  configura¬ 
tion.  Let  the  corresponding  value  of  (and,  hence,  also  that  of  o>)  be  to0°.  It  then 
follows  from  (3.19),  and  (3.2),  that  the  vorticity  vector  to  at  time  t  can  be  written 


=7fh°+ ,/*’*)■ 


The  d'Alembtrt-Euler  condition.  The  condition  that 

s-O 


(3.20) 


(3.21) 


1  *  Truesdcll  [1,  Sect.  85]  obtained  (3.20)  by  another  method  and  called  it  the  “basic 
vorticity  formula".  He  interpreted  it  as  representing  the  convection  and  diffusion  of 
vorticity.  through  the  first  and  second  terms,  respectively,  on  its  right-hand  side.  In  this 
connection,  see  also  [2]. 
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for  each  particle  and  for  all  t,  is  called  the  d'Alembert- Euler  condition.11  In  view 
of  (3.3),  (2.14),  (2.9)31,  and  (3.14)*,  each  of  the  following  conditions  is  equivalent 
to  (3.21): 

=  0,  w  =  0,  Aikw  =  0,  A,kw°  =  0.  (3.22) 

It  is  then  obvious  from  (3.19),  that  the  d’Alembert-Euler  condition  is  met  if  and 
only  if  the  time-dependent  vector  <u°  retains  the  value  co0°  which  it  had  at  time  t0, 

i.e., 

tu°  =  o>0°,  (3.23) 

and  equivalently  if  and  only  if 

W°  =  W0\  (3.24) 

i.e.,  the  value  of  W°  at  t  =  t0.  Furthermore,  it  follows  from  (3.2),  and  (3.23) 
that  a  necessary  and  sufficient  condition  for  (3.21)  to  hold  is  that 

Joj  =  Foj03.  (3.25) 

Equation  (3.25)  is  Cauchy's  vorticity  formula.  It  is  essential  to  observe  the  following 
difference  between  (3.25)  and  (3.2),:  in  the  former  a>0°  is  a  time-independent 
vector,  whereas  in  the  latter  to0  is,  in  general,  time-dependent. 

An  important  well-known  result  ( The  Lagrange-Cauchy  theorem )  can  be  im¬ 
mediately  deduced  from  (3.25):  Suppose  that  the  d'Alembert- Euler  condition  is 
satisfied.  If  a  motion  is  irrotationai  at  any  one  time,  then  it  is  irrotational  for  all 
rime.  For.  if  io(t  *)  =  0  for  some  given  instant  t*,  then,  by  virtue  of  (3.25),  a*o°  =  0. 
Hence,  again  invoking  (3.25),  we  find  that  <o(f)  =  0  for  all  t. 


4.  Circulation:  Eulsrian  a ad  Lagrangian  Forms 

Consider  a  simple  closed  curve  ct,  which  we  shall  refer  to  as  a  circuit,  in  the 
present  configuration  n,  and  let  denote  its  inverse  image  (also  a  circuit)  in  the 
reference  configuration  *°.  The  circulation  around  ‘i  has  the  representations 

r  =  f(«,  t)  =  J  vdx  = 

where  (3.11),  has  been  used.  Applying  Stokes’s  theorem  to  (4. 1)2,3.  aQd  recalling 
(2.4a)j  and  (3.15)*,  we  obtain 

at  nda  =  jta°  •  NdA,  (4.2) 

where  S'  is  any  surface  bounded  by  <4,  n  is  a  unit  normal  to  S'  (chosen  to  be  com¬ 
patible  with  the  orientation  of  •#),  -S'0  is  the  inverse  image  of  S'  in  *°,  and  N  is  a 
unit  normal  to  S'0.  It  is  clear  from  (4.2),  that  a  motion  is  irrotational  at  time  t 
if  and  only  if  f  =  0  for  every  reducible11  circuit  V  that  can  be  drawn  in  the  present 


11  See  Truesdeu.  [t,  p.  87]. 

11  i.e.,  one  that  can  be  continuously  shrunk  to  a  point  of  the  present  configuration 
while  remaining  throughout  in  this  configuration. 
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configuration  of  the  continuum  (. Kelvins  Kinematical  Theorem).  A  similar  statement 
involving  “if0  follows  from  (4.2)2. 

If  we  now  fix  attention  on  the  particles  that  compose  then  clearly  a  material 
curve  is  described  as  time  progresses  and  its  image  in  m  is  <€.  The  rate  of  change 
of  r  for  this  fixed  set  of  particles  is  given  by 

df 

r=—(V°,t)  =  Jodx=Ja°dX,  (4.3) 

where  (3.1)2  has  been  used  in  deriving  the  last  step. 

Again,  applying  Stokes’s  theorem  to  the  integrals  in  (4.3)3  4  and  recalling  (2  9)2, 
(2.14).  (3.17),,  and  (3.19),.  we  find  that 

f  =  j  x  n  da  =  j  at  n  da  =  /  s®  ■  N  dA  =  /  o»°  •  N  dA.  (4.4) 
y  y  y-  y- 

A  motion  is  circulation-preserving  if  and  only  if  T  =  0,  for  all  r,  and  for  every 
reducible  circuit  V °  that  can  be  drawn  in  the  reference  configuration  of  the  con¬ 
tinuum.  It  is  obvious  from  (4.4),  and  (3.21)  that  a  motion  is  circulation-preserving 
if  and  only  if  the  d’Alembert-Euler  condition  is  satisfied  [I,  Sect.  46]. 

For  a  a=  0,  a  vortex-line  in  the  present  configuration  is  a  line  which  is  every¬ 
where  tangent  to  the  vector  cu;  it  may  be  called  an  cu-line.  It  passes  through 
a  lefinite  set  of  particles14  in  the  configuration  x  at  time  t.  Correspondingly,  since 
dX  -  F  '1  dx,  it  is  dear  from  (3.2),  that  the  line  joining  this  same  set  of  panicles 
in  x°  is  tangent  to  the  vonicity  vector  cu®,  so  that  it  is  an  a>°-line. 

In  the  present  configuration,  draw  any  reducible  circuit  V  whose  tangent  at 
every  point  r  .  ot  r  *ai  “  to  uie  vorticity  vector  at  that  point.  Thetu-lmes  passing 
through  <<£  then  form  a  vortex-tube1*  in  *;  let  us  call  it  an  co-tube.  Transport¬ 
ing  tlii j.  tube  back  into  w\  we  obtain  an  co®-tube  there  (Figure  1). 

Let  oj(>  0)  and  <u°(>  0)  be  the  magnitudes  of  at  and  o>®,  respectively. 


and  write 

o»  =  cum, 

(4.5) 

Recall  the  formulae 

Am  =  Fm°,  AJ  =  m°  •  Cm®, 

(46) 

where  A  >  0  is  the  stretch  of  a  line-element  lying  along  the  unit  vector  m°  in  w3 
and  where  use  has  been  made  of  (2.1)5.l#  It  follows  from  (3.2),,  (4.5), j2,  and  (4.6), 


14  Of  course,  a  moment  later,  the  vortex-line  and  the  line  through  the  particles  will  in 
general  no  longer  be  identical  the  vortex-lines  are  not,  in  general,  material  lines. 
Nevertheless,  the  above  construction  is  useful.  Conditions  under  which  the  vortex-lines 
are  material  will  be  established  in  Sect.  5. 

1 9  Several  slightly  different  definitions  of  a  vortex-tube  have  been  given  in  the  litera¬ 
ture.  The  dreuit  is  sometimes  taken  to  be  infinitely  small  (as  in  Helmholtz’s  [11]  original 
definition  and  also  in  Lamb  [6]).  Redudbility  is  often  only  implteit,  and  occasionally  not 
made  part  of  the  definition  at  all  (see  Meyer  [18]).  Here,  we  have  adopted  the  definition 
given  in  Batchelor  [9]. 

1 *  The  Jcakr  function  A  is  a  function  of  the  variables  (X,  t,  m°)  in  its  Lagrangian 
representation  and  a  function  of  (*,  r,  m>  in  its  Eulerian  representation. 
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fig-  1  -  Lagrangian  and  Eulerian  vortex-tubes 


that 


0) 


(4.7) 


As  indicated  in  Figure  1,  let  y°,  &  be  tube  sections  bounded  by  circuits  <t°. 

respectively,  and  let  N,  n  be  unit  normals  vector  fields  to  these  surfaces  17. 
Then,  remembering  that 

ytt  =  F*S,  (4.8) 


where  y  (=  da  id  A)  is  the  area-stretch,  and  again  using  (3.2),.  we  have 


to  -  n  —  — a»°  •  N. 

y 


(4.9) 


Let  a  volume  dement  dV  in  be  constructed  with  its  base  dA  on  the  section 
y°.  and  with  its  sides  of  length  dL  composed  of  a*°-lines.  With  these  notations, 
dV  =  dL  dA  m°  N,  and  its  image  in  n  is  the  volume  element  dv  =  dl  da  m  ■  n. 
Since  dv  =  J dV.  it  is  then  clear  that 


J 

m  ■  n  =  —  m°  •  N. 


(4.10) 


If  is  a  normal  cross-section  (N  —  m°)  of  the  to^-tube  element,  then  m  ■  n 
=  Jiky. 


Using  (2.5)*  and  applying  the  divergence  theorem  to  a  region  in  the  present  con¬ 
figuration  bounded  by  an  <o-tube  and  two  cross-secdons  if  and  !/,  we  obtain 

Helmholtz's  First  Vortidty  Theorem. 

Jot  •  n  da  —  J  to  •  n  da,  (4  11) 


The  typeface  in  the  figures  differs  somewhat  from  that  in  the  text. 
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i.e.,  the  strength  of  an  co-tube  is  the  same  at  all  cross-section*  of  the  tube  at 
time  t  [6]. 

Equivalently,  in  view  of  (4.2)!,  we  have  Kelvin’s  form  of  this  theorem1*: 

o  =  m,  /),  (4.i2) 

where  the  circuit  V  is  the  boundary  of  the  section  y.  It  is  clear  that  the  theorem 
also  holds  for  coJ-tubes: 

f  tu°  ■  N dA  =  J  tu°  ■  N  dA,  r)  =  f(Va,  t).  (4.13) 

J'  y. 

[n  addition,  it  should  be  observed  that  since  the  vector  fields  %  and  also  are 
divergence -free,  Helmholtz’s  First  Theorem  also  holds  when  expressed  with  reference 
to  tubes  of  these  vector  fields. 

As  was  mentioned  previously,  vortex-tubes  are  not,  in  general,  material  tubes. 
Consider,  for  instance,  the  set  of  particles  that  lie  on  an  to-tube  at  time  t  (Figure  2). 
Then,  at  some  later  time  t\  the  particles  which  were  on  the  curves  V(t)  and  'if(f) 
will  now  lie  on  curves  ^(t  )  and  V(t'),  and  the  particles  on  the  tube  MNPQ  will  be 
deformed  into  the  tube  M’N'P'Q’.  However,  at  time  t\  the  vortex-tube  through 
<€(t')  could  be  M'N’P"Q"  rather  than  M'N'P'Q’,  and  the  vortex-tube  through 

£(/*)  could  be  P'Q'M'N". 


Fig.  2.  The  panicles  which  lie  on  a  vortex-tube  MNPQ  at  time  t  do  not  in  general  i:s 
on  a  vortex  tube  at  a  subsequent  time  t\ 


We  now  proceed  to  examine  the  rate  at  wh;ch  vortex-lines  change  for  a  given 
particle.  To  this  end.  we  first  take  the  material  derivative  of  (4.6),  and  employ 
(2.3),  and  (2.5),  to  obtain 


Am  4-  Am  =  A  Dm  -r  Fm°. 


(4.14) 


1  *  See  Kelvin  (12).  For  historical  comments,  see  Truesdell  [1], 
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Since  m  is  a  unit  vector,  it  follows  from  (4.14)  that 
X  m°  Cm° 

m  =  |d  —  |m  •  Dm  +  m°  •  Cm0j  /j  m  -j-  y  Fm°, 

where  use  has  been  made  of  (2. 1)3  and  (4.6),  2. 

By  virtue  of  (4.7), 


a>  a>°  X  J 


0)  CD 


Also,  in  view  of  (4.5)„  (2.14),  and  (2.13), 

t 

%  =  com  +  com 

and  hence,  by  (3.3), 

%°  =  {m}. 

It  then  follows  from  (4.18)  that 

(O  I  x  „  /  Oi  (l)°\ 

—  -t-1  m  =  ffl°  - ;)  m  , 

CO  \<U  CO  / 

where  (3.19),,  (4.5)j,  (4.6),,  and  (4.7)  have  been  utilized. 


(4.15) 


(4.16) 


(4.17) 

(4.18) 


(4-19) 


5.  Material  vortex-lines 

A  vector-field  u  is  material  if  and  only  if 

u(*,  /)  =  Fu°(X),  (5.1) 

where  u(x,  t)  is  the  value  of  u  at  a  point  *  in  the  present  configuration,  and  u° 
is  the  value  of  u  at  X  in  the  reference  configuration,  the  vectors  x  and  X  being 
related  through  the  motion  %  of  the  continuum.  We  emphasize  that  u°  in  (5.1)  is 
independent  of  time. 

It  is  obvious  from  (5.1)  and  (2.3),  that  if  u  is  material,  then 

u  ■  Lu.  (52) 

Let  us  now  consider  (2.12);  it  is  clear  that  in  general  the  vorticity  field  is  not  ma¬ 
terial.  (But,  in  the  special  case  in  which  the  motion  is  isochoric  (/  =  1)  and  the 
d’Alembert-Euler  condition  (3.21)  is  met,  it  can  be  seen  from  (3.25)  that  the  vor¬ 
ticity  field  is  material.)  Of  wider  significance  is  the  case  in  which  only  the  vortex¬ 
lines  are  material. 
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Theorem  I.  Suppose  a*  =*=  0.  Each  of  the  following  conditions  (a),  (b),  (c),  (d)  is 
necessary  and  sufficient  for  the  oj- lines  to  be  material: 


(a) 

m°  =  0, 

(c) 

o> 

0  O 

a  =  — to, 

CO 

(b) 

i  ( io  u>°\ 

(d) 

1  CO0 

m  =  -  1 - d  m, 

a  =  oj  =  — 

\CO  co  / 

co 

Proof.  First  observe  that  the  four  conditions  in  (S.3)  are  all  equivalent  to  one  an* 
other.  Thus,  by  (4.19),  (4.6),,  and  (4.7),  condition  (S.3a)  is  equivalent  to  condition 
(5.3b).  Utilizing  (4.5),  and  (3.19),,  we  see  that  conditions  (5.3a)  and  (5.3c)  are 
equivalent  to  each  other.  Finally,  by  using  the  transformations  (3.2),  and  (3.3) 
and  their  inverses,  we  can  establish  the  equivalency  of  (5.3c)  and  (5.3d). 

(i)  Necessity.  By  the  definition  (5.1),  the  vector  m°  in  (4.6),  must  be  independent 
of  time.  Therefore,  condition  (5.3a)  holds.  Hence  also  do  (5.3b),  (5.3c),  and 
(5.3  d). 

(ii)  Sufficiency.  If  condition  (5.3a)  holds,  then 

m°  =  m0°,  (5.4) 

a  time-independent  unit  vector.  Substituting  (5.4)  into  (4.6),,  we  obtain 

Xm  —  Fm0°,  (5.5) 

which  implies  that  the  oi-lines  are  material.  Since  we  have  already  shown  that  each 
of  the  conditions  (5.3  b),  (5.3c),  and  (5.3d)  is  equivalent  to  (5.3a),  they  also  must 
individually  imply  (5.5). 

Corollary  I  (Helmholtz’s  Second  Vortidty  Theorem).  If  the  d' Alembert- Euler  con¬ 
dition  is  satisfied ,  the  vortex-lines  are  material  lines. 

An  equivalent  statement  is:  If  the  motion  is  circulation-preserving,  the  vortex¬ 
lines  are  material  lines.  To  prove  the  corollary,  simply  observe  that  (3.22),  is 
equivalent  to  the  d’Alembert-Euler  condition,  and  consequently,  in  view  of  (3.19), 
and  (4.5),,  m°  —  0.  Therefore,  condition  (5.3  a)  is  satisfied,  and  the  result  follows 
from  the  theorem. 

It  is  of  interest  to  note  that  in  the  present  case  (i.e.,  when  a0  =  0),  (5.3b) 
yields 

t  OJ 

m  - - m.  (5.6) 

a» 


Corollary  O  (Helmholtz’s  Third  Vortidty  Theorem).  If  the  d'Alembert- Euler  con¬ 
dition  is  satisfied,  a  time-independent  strength  can  be  associated  with  each  vortex-tube. 

Equivalently,  we  may  replace  the  hypothesis  of  this  corollary  with  the  hypo¬ 
thesis  that  the  motion  be  circulation-preserving.  The  proof  is  straightforward: 
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Helmholtz’s  Second  Theorem  implies  that  the  vortex-tubes  are  material,  and  (4.4), 
now  tells  us  that  f  =  0.  Therefore,  in  view  of  (4.1)4, 

r=f(V°).  (5.7) 

Hence,  recalling  Helmholtz’s  First  Theorem,  we  see  that  a  time-independent 
strength  can  be  associated  with  each  vortex-tube. 

The  upshot  of  the  three  Helmholtz  Theorems  is  that  if  the  d'Alembert-Euler 
condition  is  satisfied,  then  each  vortex-tube  moves  as  a  material  surface,  and  with 
it  is  associated  a  scalar  measure  r  which  varies  neither  along  the  tube  at  a  given 
value  of  t,  nor  as  a  function  of  t  for  a  given  tube.  Thus,  once  the  d’Alembert- 
Euler  condition  is  satisfied,  r  is  a  constant  for  each  vortex-tube. 

Corollary  m.  If  the  ut-lines  are  material  and  if 

oi°  =  0,  (5.8) 

then  the  d'Alembert-Euler  condition  is  satisfied. 

This  follows  immediately  from  (5.3a),  (4.5)2,  (3.19),,  and  (3.22),. 

Corollary  IV.  In  a  circulation-preserving  flow,  no  stretching  of  vortex-lines  occurs 
if  and  only  if  Jui  is  time-independent.19 

In  a  circulation-preserving  flow,  <*  =  0,  and  hence  by  (3.23),  <o°  =  <o0o. 
Also,  the  vortex-lines  are  material  by  Corollary  I.  Setting  k  =  1  in  (4.7),  we  then 
see  that  J<o  =  o>0°,  which  is  time-independent. 

Theorem  O.  Suppose  <0*0.  Each  of  the  following  conditions  (a),  (b),  (c),  (d) 
is  necessary  and  sufficient  for  the  ut-lines  to  be  material : 


(a) 

m  x  m°  =  0, 

(b) 

<o°xco°  =  0, 

(c) 

to  x  *  =  0, 

(d) 

mx  m  =  0. 

Proof.  It  follows  from  (3.2),,  (3.3),  and  (3.  l)i,2  that 

cox*  = -j(F~l)r(co°x»°).  (5.10) 

Therefore,  cox*  =0  if  and  only  if  <o°x<*°  =  0.  By  virtue  of  (3.19),  and  (4.5)2, 
the  latter  can  occur  if  and  only  if  m°  x  m°  =  0,  which  can  occur  if,  and  (since 
m°  is  a  unit  vector)  only  if,  m°  =  0.  Also,  in  view  of  (4.17)  and (4.5),,  <ox»  =  0 

if  and  only  if  m  x  m  =*  0.  Consequently,  each  of  the  conditions  in  (5.9)  is  equi¬ 
valent  to  (5.3a),  which  is  a  necessary  and  sufficient  condition  for  the  co-lines  to  be 
material.20 

■’  [I.  p.  180]. 

i0  The  condition  (5.9c)  is  given  by  Truesdell  [1,  Sect.  45],  who  employs  a  different 
argument. 
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